Monte Carlo computer simulations and electron microscopy of colloidal cluster 
formation via emulsion droplet evaporation 
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We consider a theoretical model for a binary mixture of colloidal particles and spherical emulsion 
droplets. The hard sphere colloids interact via additional short-ranged attraction and long-ranged 
repulsion. The droplet-colloid interaction is an attractive well at the droplet surface, which induces 
the Pickering effect. The droplet-droplet interaction is a hard-core interaction. The droplets shrink 
in time, which models the evaporation of the dispersed (oil) phase, and we use Monte Carlo sim- 
ulations for the dynamics. In the experiments, polystyrene particles were assembled using toluene 
droplets as templates. The arrangement of the particles on the surface of the droplets was analyzed 
with cryogenic field emission scanning electron microscopy. Before evaporation of the oil, the particle 
distribution on the droplet surface was found to be disordered in experiments, and the simulations 
reproduce this effect. After complete evaporation, ordered colloidal clusters are formed that are sta- 
ble against thermal fluctuations. Both in the simulations and with held emission scanning electron 
microscopy, we find stable packings that range from doublets, triplets, and tetrahedra to complex 
polyhedra of colloids. The simulated cluster structures and size distribution agree well with the 
experimental results. We also simulate hierarchical assembly in a mixture of tetrahedral clusters 
and droplets, and find supercluster structures with morphologies that are more complex than those 
of clusters of single particles. 

PACS numbers: 61.46.Bc,61.20. Ja,68.37.-d 



I. INTRODUCTION 

Systems of colloidal particles that interact with pair 
potentials that feature short-ranged attraction have been 
extensively studied in the literature. The typical bulk 
phase diagram shows a characteristic gas-liquid phase 
transition that can become metastable with respect to 
freezing for sufficiently short range of attraction [THll- 
While the addition of a short-ranged soft repulsion is 
known to change the phase diagram Q only quantita- 
tively, a repulsion with a range larger than that of the at- 
traction can have a dramatic effect on the topology of the 
phase diagram 0|. Moreover, both in experiments and 
in dynamical simulations the equilibrium phase diagram 
for systems with competing interactions is often over- 
shadowed by non-equilibrium phenomena such as vitrifi- 
cation, gelation and cluster formation (sj-fl^. Typically, 
pair potentials with a short-ranged attraction show spon- 
taneous clustering of particles at sufficiently low tempera- 
tures. The geometric structure of clusters of particles has 
been studied theoretically, e.g. for the Lennard- Jones po- 
tential [l^, , hard spheres [l^ , and hard spheres that 
additionally interact with a short-ranged attraction [l6j . 
The structures obtained using the Morse potential were 
analyzed by Doye, Wales, and Berry and by Taffs 
et al. [3 and compared to those from the Asakura- 
Oosawa potential H^]. Stable clusters of colloids 
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are interesting because they can be viewed as colloidal 
molecules [21,, ,221 that can potentially be used as build- 
ing blocks for the fabrication of novel materials. The size 
and geometric structure of colloidal clusters, however, are 
not easily controllable in experiments. Several different 
methods that offer control of the clustering process have 
been proposed. Jiang and Granick prepared clusters 
using Janus colloidal particles, i.e. spherical particles 
that possess oppositely charged hemispheres. Cluster of 
particles with larger numbers of patches have also been 
studied (23 - [26| . Erb et al. (23| succeeded in preparing 
clusters of magnetic particles. 

One particularly promising approach, based on the 
evaporation of the dispersed phase in an emulsion, 
was developed by Velev, Furusawa, and Nagayama 
(29I . [soj . Here colloidal particles adsorb at the interface 
between dispersed and continuous phase in order to min- 
imize the interfacial free energy (Pickering effect) (3l| . 
During evaporation the particles are pushed together by 
capillary forces and subsequently held together by van 
der Waals interactions. Manoharan, Elsesser, and Pine 
(3^ prepared micron-sized clusters using this technique 
with polystyrene microspheres that were 844 nm in diam- 
eter. The authors found clusters of particles with pack- 
ings that minimize the second moment of the mass dis- 
tribution. The emulsion method is versatile and was sub- 
sequently used to obtain clusters of particles 220 nm in 
diameter [s^ , of patchy particles [SJ, , and of bidis- 
perse colloids [s^. Shear was used by Zerrouki et al. (37| 
to produce monodispersed droplets. A similar technique 
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is based on aerosol droplets [38| instead of oil droplets. 

Similar in spirit to the emulsion evaporation technique, 
an alternative miniemulsion technique (sol . l40j was re- 
cently developed. Here, a miniemulsion is prepared from 
a standard emulsion by ultrasonication. The sound waves 
produce an emulsion of small droplets in a process of fis- 
sion and fusion [65]. The average size of the droplets can 
be tuned in the range of 360 nm to 1800 nm kll. Small 
colloidal particles can then be used in solution with the 
small monodispersed droplets of the miniemulsion to ob- 
tain clusters |40| that have diameters much smaller than 
1 fim. Although these clusters can consist of many con- 
stituent colloidal spheres, they remain small enough to 
reside in the colloidal domain [42|, i.e. Brownian mo- 
tion thermalizes such systems. Therefore hierarchical 
self-assembly comes within reach. 

In contrast to the large body of experimental work and 
the closely related theoretical efforts to understand the 
resulting cluster structures and their symmetries, little 
theoretical work has been done to describe the process of 
cluster formation. Roman, Schmidt, and Lowen pro- 
posed a model for a dispersion of hard spheres and emul- 
sion droplets, but these authors did not investigate clus- 
ter formation. Lauga and Brenner (4^ modeled and sim- 
ulated the evaporation-driven assembly of colloidal par- 
ticles. They considered individual droplets with varying 
numbers of adsorbed particles and calculated the (non- 
spherical, in general) shape of the oil-water interface by 
the requirement of minimal surface free energy. They 
considered different values of the contact angle and re- 
ported good agreement with experimental findings. Very 
recently, Mani et al. [4^ studied the stability of larger 
coUoidosome-like shells of particles, albeit without mod- 
eling the assembly process. 

In this paper we present a basic model to describe 
the process of cluster assembly through emulsion droplet 
evaporation. We use Monte Carlo computer simulations 
to study the cluster formation of colloids with competing 
short-ranged attraction and long-ranged repulsion inter- 
actions. Differently from Lauga and Brenner [ZJl we also 
simulate the dynamical capture of the colloidal particles 
onto droplet surfaces and studied not only the cluster 
structures, but also analyzes the histograms of the clus- 
ter size distribution. We complement these calculations 
with experiments of polystyrene particles in an oil-in- 
water emulsion. We use cryogenic field emission scan- 
ning electron microscopy (cryo-FESEM) and field emis- 
sion scanning electron microscopy (FESEM) to investi- 
gate the intermediate and final stages of cluster forma- 
tion, respectively. The cryo-FESEM micrographs show 
the distribution of small particles on a droplet surface. 
We find the simulation results to be in good agreement 
with the experimental results for both intermediate and 
final cluster structures. We also find good quantitative 
agreement between experimental and simulation results 
for the cluster size distribution. Having demonstrated 
the validity of our model, we study the possibility of hi- 
erarchical self-assembly, by carrying out simulations of 



a mixture of thermal tetrahedral clusters and emulsion 
droplets. We obtain clusters of clusters (superclusters) 
with structures that differ from the clusters made of sin- 
gle particles. 

The paper is organized as follows. In Sec. Ill Al we give 
the details of the pair interactions. In Sec. Ill B l and IIIC] 
we present simulation and experimental details, respec- 
tively. In Sec. nil Al we describe the results for the dy- 
namics of cluster formation. In Sec. IIII Bl we show the 
results for the cluster structures and for the histograms of 
the size distributions. In Sec. IIII Cl we present the results 
for the superclusters. Final remarks and conclusions are 
given in Sec. IIVI 



II. MODEL AND METHODS 

A. Definition of the pair interactions 

We study a binary mixture of Nc colloidal particles 
with hard-sphere diameter CTc and Nd droplets of diam- 
eter Ud- The total interaction energy U is the sum of 
colloid-colloid, droplet-droplet and colloid-droplet inter- 
actions, 

iJ = ^ (/),,(|r, - r,|) + ^ 0,rf(|R, - R,|) 

i<j i<j 

where r.^ is the center-of-mass position of colloid i, Rj is 
the center-of-mass position of droplet j, (j)cc is the colloid- 
colloid pair interaction, (f>cd is the colloid-droplet pair in- 
teraction, and 4'dd is the droplet-droplet pair interaction. 

We consider two different types of colloid-colloid in- 
teractions. The first is the sum of the short-ranged at- 
tractive Asakura-Oosawa potential UAoif) ^^nd the long- 
ranged repulsive Yukawa potential UY{r), i.e. 

, / X / oo r < ac 

9cc[r) - I jj^^^^ ^ UAoir) otherwise, 

with /? = l/fc^T, where ks is the Boltzmann constant 
and T the temperature. Here the Yukawa potential is 
defined by 



^-K(r-(Tc) 



(3) 



and the Asakura-Oosawa potential [19|, |20| is 



^Ao(^) - otherwise. 



with /(r) 



' jyr) = - 2[i+q)a, + 2(i+g)^ag J ' ^^"^ parameter 
controls the strength of attraction, while ey controls 
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Figure 1. Pair interactions for the binary mixture of col- 
loids and droplets, a) Comparison between the Asakura- 
Oosawa- Yukawa (AOY) and the Square-Well- Yukawa (SWY) 
potentials with q — 0.1, nac=10, /3ey=24.8, /Seao ~ 2.1, 
P^sw = 9- These parameters were chosen in order to obtain 
|min(/?(/!)cc)| = max(/3(/!>cc) = 9. b) Colloid-droplet potential 
4'cdHo''^ scaled by the droplet-solvent interfacial tension 7 for 
o'd/o'c = 3, 2, 1.1, and 1 (from right to left). 



the strength of the repulsion. The range of the interac- 
tions are controlled by the parameter q for the AO poten- 
tial and by k for the Yukawa interaction. We refer to the 
combined interaction ^ as the Asakur a- Oosawa- Yukawa 
(AOY) potential. Although both the AO and Yukawa 
potentials have very specific physical interpretations, we 
use them in this paper merely as generic models for a 
steep short-range attraction and a long-range repulsion. 



respectively. The AO potential describes the depletion 
attraction between colloidal particles due to the presence 
of non-adsorbing polymers with radius of gyration (Jc9/2. 
In the limit of small q Eq.Q is exact and f{r) is related 
to the free volume gained by the polymers when two col- 
loids are close to each other The Yukawa potential 
describes the interaction between two charged particles 
screened by a medium with inverse Debye length k. As 
shown in Fig. [T^), for a typical set of parameters (justi- 
fied below) that we used in the simulations, the potential 
4>cc{t) defined by Eq. ([2]) has a maximum at distance 
^poak- In order to investigate the influence of the shape 
of the attractive part of the potential we also consider a 
modified version of ([2]) using a square well potential for 
distances smaller than rpcak- This is 



Urir) 



r < (Tc 

(7r. <r < Tpoak 



UAoif) otherwise. 



(5) 



which we refer to as the Square- Well- Yukawa (SWY) po- 
tential. It is shown in Fig. [T] with a dashed line. This 
potential is very similar to that used by Mani et al. 
|45| in their pioneering study of the stability of colloidal 
shells. The parameter space of the interaction is ar- 
bitrarily restricted to potentials with the shape shown 
in Fig. [U i.e. with max(0cc) = |min((/)cc)|, and hence 
0cc(?'poak) = ^0cc(o'c)- The height of the repulsive bar- 
rier, (/)cc(?'poak) — 0cc(oo), IS half the depth of the attrac- 
tive well, (/>cc(fpoak) - </'cc(o-c)- 

The droplet-droplet interaction is taken to be hard- 
core repulsion 



(l>dd{r) 



00 r < ad + Uc 
otherwise, 



(6) 



with an effective hard-core diameter Ud + Cc that is larger 
than the bare droplet diameter Ud- Using the effective di- 
ameter ensures that the surface-surface distance between 
any two droplets is always larger than one colloid diam- 
eter. In this way two droplets can never bind together 
due to a shared colloid. 

The colloid- droplet interaction is aimed at mode ling 
the Pickering effect. The loss of interfacial energy (46| 
when a particle is trapped at the surface of the droplet 
is 76". with S the droplet surface that is covered by the 
colloid, and 7 the droplet-solvent interfacial tension. The 
relation is valid when the interfacial tension between the 
colloid and the droplet is the same as that between the 
colloid and the solvent. The surface S has different ex- 
pressions depending on the size of the droplet and the 
colloid-droplet separation. If the diameter of the droplets 
is larger than the diameter of the colloidal particles, i.e. 
for (Trf > (Tc, the colloid-droplet energy is 



otherwise, 



0'd + (^a 



(7) 



with/i = (crc/2-<Td/2-hr)((Tc/2+f7d/2-r)/(2r) the height 
of the spherical cap that results from the colloid-droplet 
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intersection. On the other hand, when the diameter of 
the droplets is smaller than the diameter of the colloidal 
particles, i.e. ad < Oc, we assume 

</.cd(0 = <^ --iirodh 2^ < r < (8) 
[ otherwise. 

Within this model we neglect the influence of the par- 
ticle on the oil-water interfacial curvature. The contact 
angle changes upon changing the position of the particle 
with respect to the oil-water interface, while in reality 
the contact angle remains constant and the curvature of 
the oil- water interface changes, such that the droplet be- 
comes non-spherical, which is beyond our model. 



B. Simulation method 

In our model we neglect the coalescence of the droplets 
and the hydrodynamic interactions due to the solvent. 
Even with these simplifications, the relevant time scale 
is not easily reachable in standard Brownian Dynam- 
ics (BD) simulations. Hence, we carry out Metropolis 
Monte Carlo (MC) simulations, that for small displace- 
ment steps can reproduce the correct dynamics |47| and 
are much more efficient to run than BD simulations. The 
evolution of the system is therefore described by the num- 
ber of MC sweeps per particle. For the colloidal parti- 
cles [l^l the MC maximum trial displacement dc is set to 
dc = O.OIctc- 

We define a bond between two colloidal particles when 
their distance is smaller than r-poak- A cluster is de- 
fined as a set of colloidal particles connected by a net- 
work of bonds. The translational diffusion of the clusters 
is achieved by translational MC cluster moves j4^ with 
maximum linear displacement = dd \/N with TV the 
number of particles in the cluster. This approximates 
the hydrodynamic slowing down of a spherical cluster 
that satisfies the Stokes-Einstein equation for the diffu- 
sion constant D = ^^^J^ , with rj the viscosity of the 
solvent and ads the diameter of the sphere that approx- 
imate the shape of the cluster. Here, we assume that 
f cis Additionally, we mimic the cluster rotational 

diffusion via rotational MC moves, in which clusters are 
rotated around a random axis with a maximum angle 
^cis = O.OIo-c/ctcIs. In a cluster move jH, |4§|, all parti- 
cles belonging to a cluster are translated or rotated by 
the same amount. In order to satisfy the condition of 
detailed balance all cluster moves that lead to two clus- 
ters merging or a cluster and a single particle merging 
are rejected. 

The droplets move according to th e MC scheme with 
a maximum displacement dd = da\/ a^j ad- The evap- 
oration dynamics is introduced by forcing the droplets 
diameter ad to shrink at a fixed rate. The rate is chosen 
so that the droplets vanish half-way though the simula- 
tion (5 X 10^ sweeps). This leaves another 5 x 10^ sweeps 



to investigate the stability of the clusters against thermal 
fluctuations. Here the timescales are chosen for practical 
reasons. In the experiments stability can be relevant on 
the time scales of years, while the clusters experiments 
last typically tens of minutes. Hence, our simulation do 
not address the true long-time behavior of the system. 

We restrict ourselves to symmetric potentials with 
max(/30cc) = |miii(/3<?I>cc)|- In particular, we investigate 
the case max(/30cc) = 9 as shown in Fig. H] (g = 0.1, 
KCTc = 10 and Tc = 2.5 (Tc)- The height of the repulsive 
barrier is 9kBT, a value big enough to hinder sponta- 
neous clustering, while the depth of the attractive well 
is ISfc^T, so that in practice a particle cannot escape by 
thermal fluctuations. For these parameters we flnd that 
the maximum of the potential is at rpeak = 1.0845(Tc. For 
each parameter set we run eight independent MC simula- 
tions with Nc ~ 500 colloidal particles with packing frac- 
tion r/c = 0.0034, and 0.01. The droplets packing fraction 
is fixed at rjd = 0.1. Simulations have been performed for 
initial droplet sizes ad{0) = 2,4,6, and Sctc- Each sim- 
ulation consisted of 10^ MC sweeps. In every sweep all 
particles are attempted to be moved on average once. 
The droplets and the particles are initialized randomly 
in the cubic simulation box, with the constraints that all 
colloids are outside of the droplets and that the mini- 
mum distance between colloidal particles is larger than 
Tpoak- Hence, we start in a state without clusters. We 
characterize the cluster structure by the bond-number 
Uf), corresponding to the number of pairwise bonds in a 
cluster. The number of bonds is also an estimate of the 
total energy of the cluster; a higher number of bonds 
corresponds to a greater attractive energy. 



C. Experimental methods 

1. Cluster preparation 

The colloidal particles are positively charged and nar- 
rowly dispersed polystyrene spheres with 154 nm diam- 
eter. A detailed description of the preparation of the 
constituent particles and their assembly into clusters is 
given in Ref. [s^- In particular, in these experiments, the 
particles were added both via the water and via the oil 
phase. For the present studies the clusters were prepared 
in a slightly modified fashion. Briefly, 53 mg polystyrene 
particles suspended in 3 ml of toluene and 73 fA dodecane 
(to suppress Ostwald ripening) were emulsified with 27 
ml of an aqueous solution of Pluronic® F-68 (1 wt.%) us- 
ing an ultrasonic homogenizer (Sonoplus HD 3200, Ban- 
delin) . Evaporation of the dispersed toluene phase under 
reduced pressure (50 mbar, 40 °C) initiated the assembly 
of the particles into clusters. 
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2. Electron microscopy of droplet and assembly morphology 

The emulsion droplets bearing polystyrene particles at 
their surface were examined on a cryogenic field emis- 
sion scanning electron microscope (Ultra Plus, Zeiss). 
Specimen preparation was accomplished by sandwiching 
4 [A of the emulsion in between two aluminum platelets 
(3 mm X 0.5 mm, 0.15/0.15 mm. Engineering Office 
M. Wohlwend). The carrier assembly was plunged into 
a high-pressure-freezing machine (EM HPMIOO, Leica) 
and was vitrified at 2000 bar within 20 ms. This helped 
sealing the sandwich so that nucleation of ice crystals 
and specimen damage were suppressed. In a cryo prepa- 
ration chamber (EM MED020 FF, Leica) the sample was 
freeze-fractured, lightly etched for 60 s at —112 °C, and 
sputtered with platinum in an amount equivalent to a 
4 nm thick coating. The specimen was transferred by a 
cryo shuttle (EM VCT 100, Leica) to the cold stage of 
the microscope. Micrographs were recorded digitally at a 
temperature of —160 °C, with an aperture of 10 /xm and 
a voltage of 1.0 kV. The morphologies of the colloidal 
assemblies were analyzed by field emission scanning elec- 
tron microscopy (FESEM) on a Zeiss LEO 1530 Gemini 
microscope equipped with a field emission cathode oper- 
ating at 3 kV. A minute amount of the cluster suspension 
(10~^ wt.%) was placed onto a silicon wafer (CrysTec) 
and dried under ambient conditions. The specimen was 
coated with a platinum layer of 1.3 nm thickness using 
a sputter coater (Cressington 208HR) to make the spec- 
imen conductive. 



III. RESULTS 

A. Dynamics of cluster formation 

Figure [5] shows simulation snapshots at four differ- 
ent stages of the simulation. The initial configuration 
of the simulation (Fig. (2^-) is a binary mixture of non- 
overlapping spheres. The large (pink) spheres represent 
the droplets, while the small (blue) spheres represent 
single colloidal particles. After 3.6 x 10'^ MC sweeps 
(Fig. [2h) particles are trapped at the surface of a droplet 
(red). Figure [2t) shows the configuration after 5 x 10^ 
sweeps, just after the droplets have vanished completely. 
Figure [2Ji) shows the configuration at the end of the sim- 
ulation after 10^ sweeps. All clusters that are formed 
formed due to the droplets are still present in the system 
demonstrating the stability of the clusters against ther- 
mal fiuctuations. Few doublets (green) have also formed 
spontaneously. ISee also a movie of the simulation here.l 
As a check we simulated the structure of the pure col- 
loidal fiuid without droplets. We find that the structure 
of the single component fluid is largely composed of single 
particles with only 1% of particles belonging to doublets. 
Larger clusters are not formed during the span of our sim- 
ulations. This can be explained easily by the following 
argument. The minimum distance between colloidal par- 



ticles is at the beginning of the simulation larger than 
''peak- 111 order to form a bond, the colloidal particles 
have to overcome the repulsive energy barrier (OfesT) of 
the colloid-colloid interaction. The probability to ther- 
mally overcome this barrier is very small, and in order 
to form larger clusters the particles have to be forced be- 
yond the repulsive barrier. Hence, in the restricted time 
interval that is accessible in the simulations, the colloidal 
fluid is (quasi-)stable. However, for longer times further 
clustering might occur (see discussion in the Sec. IIIB[) . 

In the binary mixture, a clustering mechanism is pro- 
vided by the shrinking droplets. The formation of small 
clusters implies that only a limited number of particles 
are bound onto the droplets. The small size of our parti- 
cle makes them ideal for the self-assembly of small clus- 
ters that are well suited as building blocks for subse- 
quent self-assembly, but prevented us to follow the assem- 
bly process in experiments in real space like Manoharan, 
Elsesser, and Pine [s^l- Therefore, we studied their distri- 
bution in the experimental emulsions with cryo-FESEM. 
The micrographs (Fig. [3]) indicate a random distribu- 
tion of the positions of the particles at the droplet sur- 
face. The void in the center presents the imprint of the 
frozen dispersed phase, i.e. a single micron-sized toluene 
droplet. The polystyrene particles left in the cavity after 
sublimation of toluene are randomly distributed at the 
former droplet surface. Because toluene is a good solvent 
for polystyrene, the particles are significantly swollen at 
the droplet interface as expressed by a larger diameter. 
This agrees with the fact that the micrograph indicates 
that the polystyrene particles prefer the dispersed toluene 
phase. The position distribution of the particles on the 
droplet surface in simulations (Fig. and inset Fig. 
is in good agreement with the experiment and indicates 
that the particles can freely diffuse on the surface of the 
droplets. 

In simulations, we consider the case that the particle 
surface has no preference of whether it wets oil or water. 
Hence the contact angle at a (planar) oil-water interface 
is 90° because of the assumption that the colloid-solvent 
interfacial tension is equal to the colloid-droplet interfa- 
cial tension - see the discussion above Eq.© . Despite 
the difference in contact angle between simulation and 
experiments we do not expect the contact angle to affect 
the final results Q. 

The liquid structure can be further characterized in 
computer simulations. We calculated the colloid-droplet 
radial distribution functions, gcd{f), and the colloid- 
colloid radial distribution functions, gcc{f), at different 
stages of the simulations. Since the droplet diameter ad 
changes continuously during the simulation the resulting 
transient structures captured by the distribution func- 
tions are not at equilibrium. Fig. U] shows gcdif) and 
gccir) at different stages. For the colloid-colloid correla- 
tion function we plot both the results of the SWY (blue 
full line) and AOY potentials (yellow dashed line). In 
particular. Fig. 5^) shows the correlation functions after 
10^ MC sweeps. The colloid-droplet radial distribution 
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Figure 2. Snapshots of the binary mixture of colloids and droplets at colloid packing fraction r]c = 0.0034 and SWY potential. 
Shown are results at different stages of the simulation, a) Initial configuration, b) after 3.6 x 10^ , c) after 5 x 10^ , d) after 10® MC 
sweeps. Droplets are shown as pink spheres, single colloidal particles are depicted as blue spheres, green is used for particles 
belonging to spontaneously formed clusters. Colloidal particles trapped at the surface of a droplet or in a droplet-induced 
cluster are shown in red. A movie of the simulation is available ihere.l 



function gcdif) has a small peak at r ~ 3.1 CTc correspond- 
ing to the instantaneous droplet radius (Jd{t)/2. The peak 
is due to colloidal particles trapped at the droplet sur- 
face. At the same time the colloid-colloid radial distri- 
bution functions gcci^) are apparently flat outside of the 
core region. In fact, the radial distribution function is 
decaying in a way consistent with the Boltzmann factor 
(7(7') = p-f^'Pcair) expected for an equilibrium low den- 
sity gas. FiguredJa) shows the correlation functions after 
3 X 10^ MC sweeps. The cross pair correlation function 
gcd{T) shows that the droplets have shrunk to a radius 
ad{t)/2 ~ 1.45 (Tc, while gcc{r) has developed structure 
at intermediate distances. These results can be explained 
by particles trapped at the surface of the droplets and 
interacting with each other via the long-ranged colloid- 
colloid repulsion. After 3.6 x 10^ MC sweeps (Fig. 21:) the 
droplet radius has become smaller than CTc- The strong 
peak at r = CTc indicates that a large number of bonds 



between colloidal particles have formed for both the AOY 
and SWY potentials. Finally, after 5 x 10^ MC sweeps 
(Fig. UJl) a dramatic change of gcd{r) is observed, which 
is due to droplets having a diameter ad{t) = 0, and dif- 
fusing freely. We show the left panel of Fig. for com- 
pleteness, but stress that it does not correspond to any 
physical situation. On the other hand, gcc{i") shows that 
strong peaks have formed also at distances larger than 
(Tc, indicating the presence of small clusters. The time 
evolution of gccir) for the two potentials is very simi- 
lar. Most notably the final configurations differ due to 
different final cluster structures. 



B. Cluster structure and size distribution 

An overview of the different cluster structures found at 
the end of the simulation runs and in the experiment is 
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a) (7rf = 6.2 (7, 



Figure 3. Colloidal particles trapped at the surface of 
an emulsion droplet obtained with a Cryo-FESEM micro- 
graph of a toluene-in-water emulsion stabilized by crosslinked 
polystyrene particles. Inset: Simulation snapshot of a single 
droplet and colloidal particles trapped at its surface obtained 
after 3 x 10^ MC sweeps. 



presented in Figs. [S] and [HI In particular, Fig.[S]shows the 
simulation structures obtained at colloid packing fraction 
ijc = 0.0034 for the SWY potential and (7^(0) = Sctc- 
We find clusters with sizes between rtc = 2 (doublets 
not shown) and ric = 9, with ric the number of particles 
belonging to a cluster. 

In particular, we find that for Uc < 7 the clusters have 
the same structures as Lennard- Jones clusters We 
find triplets for ric — 3 (Fig. [5^) , tetrahedra for ric = 4 
(Fig. [SJd), triangular dipyramids for ric = 5 (Fig. Et), 
octahedrons for ric = 6 (Fig.[5ji), and pentagonal dipyra- 
mids for ric ~ 7 (Fig. [S^). For large number of parti- 
cles we find the snub disphenoid for Uc = 8 (Fig. [Sf) 
and triaugmented triangular prism for ric ~ 9 (Fig. [S^, 
with two different orientations). In computer simulations 
clusters with particle numbers ric > 10 are obtained at 
higher packing fractions. The additional structures ob- 
tained at colloid packing fraction rjc ~ 0.01 are shown 
in Fig. [S] The square dipyramid is found for ric = 5 
(Fig. [5^) Also, for the clusters with rtc = 10 ( gyroelon- 
gated square dipyramid in Fig. \Ejp) and ric = 12 (icosa- 
hedron in Fig. [BJi) we find good agreement with exper- 
iments. The cluster with ric = 11 (icosahedron minus 
one, Fig.[6j:) was, on the other hand, not found in exper- 
iments. As this structure is identical to the icosahedron 
except for one missing particle, it can easily be missed in 
the experimental FESEM micrographs. 

These structures are also in good agreement with those 
observed in previous experiments |32| . l40j . As noted 
by Manoharan, Elsesser, and Pine |32| . clusters contain- 
ing ric = S (snub disphenoid), ric = 9 (triaugmented 
triangular prism) and ric = 10 (gyroelongated square 
dipyramid) particles are members of a set of convex poly- 
liedra |50,] . 

Larger cluster structures can be obtained not only by 
increasing the colloid packing fraction but also by in- 
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Figure 4. Colloid-droplet (left column) and colloid-colloid 
(right column) radial distribution functions, Qcdir) and gcc(r), 
respectively, as a function of the scaled distance r/(jc at col- 
loid packing fraction r]c = 0.0034. We plot gcc(r-) for the SWY 
(blue full line) and the AOY (yellow dashed line) potentials. 
Shown are results at different stages of the computer simula- 
tion: a) after 10^, b) after 3 x 10^ c) after 3.6 x 10^ d) after 
5 X lO'^ MC sweeps. 
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Figure 5. Cluster structures found in simulations (left) for 
SWY potential, (7^(0) = 8(Tc, »7c ~ 0.0034 and micrographs 
from FESEM (right). The scalebars indicate 200nm. The 
wireframe in the simulation structures connects the particles 
centers in order to visualize the geometric arrangement. 



a) nc=5 Square dipyramid 



c) nc=ll Icosahiedron minus one 




b) nc=10 Gyroelongate square dipyramid d) nc=12 Icosahedron 




Figure 6. Same as Fig. O but for rjc = 0.01. Only additional 
structures not found at rjc = 0.0034 are shown. 



creasing the starting droplet size, as demonstrated by 
the histograms of the number of clusters iV„^ , with Uc 
the number of colloidal particles forming the cluster. Fig- 
ure [7] shows the cluster distribution for varying starting 
droplet sizes 2 < ad{0)/ac < 8. From Fig. Eli) to Fig.EJi) 
the distribution becomes broader for larger droplet diam- 
eters (Trf(O), while at the same time the yield of smaller 
clusters decreases. The presence of a greater number of 
large clusters at larger droplet diameters can be explained 
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Figure 7. Distribution of the number of clusters, N„^, as a 
function of the number of particles in the cluster, iic, found 
in computer simulations (SWY, rjc = 0.0034). For clarity, 
the Tic — 1 bar is omitted. Shown are result for different 
initial droplet sizes: a) (7^(0) /(Jc = 2. b) crd{0)/ac = 4. c) 
crd(0)/crc = 6. d) ad{0)/ac = 8. 



by the larger surface available to capture colloidal parti- 
cles in the initial stages of the simulation. Likewise, the 
probability to capture a small number of colloidal parti- 
cles decreases with increasing droplet surface, which leads 
to a decrease in the yield of small clusters. 

The comparison with experiments can be made more 
quantitative by comparing the experimentally (40| mea- 
sured weight fraction of the clusters with the total num- 
ber of particles, UcNn^ , belonging to a cluster with ric col- 
loids. In order to compare the two quantities we normal- 
ize the experimental results with the weight fraction of 
single particles and the simulation results by the number 
of single particles A^i. In the first experiments by Wag- 
ner et al, the particles were dispersed in the oil phase. In 
subsequent work Wagner et al. [ioj , experimentally com- 
pared this situation with that of adding particles via the 
water phase. They found that the same cluster structures 
result and that the cluster size distributions are also sim- 
ilar. In the present paper, for consistency, in Fig. 8 we 
compare the simulation results to experiments where the 
particles were added via the water phase. We keep, as 
an additional data set, the size distributions that were 
obtained by adding particles via the oil phase. Figure [8] 
shows the two experimental results together with the sim- 
ulation results for the SWY potential, (Jd{Q)/(7c = 4, and 
77c = 0.0034 corresponding to a percentage of particles 
per oil of 3.4%. As the degree of polydispersity of our 
emulsions is low [i^, we do not expect it to have a sig- 
nificant effect on the experimental size distribution, and 
the comparison to the (monodispersed) simulation results 
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increasing number of constituent particles the number of 
isomers increases. Fig. ^jp) shows the stacked histogram 
of the number of clusters with a specific bond-number 
for the AOY potential (crd(O) = 8ctc, ??c = 0.0034). Strik- 
ingly, the AOY potential produces a larger number of 
different isomers than the SWY potential. In particu- 
lar, the AOY potential has isomers with smaller bond- 
numbers compared to the SWY potential. We interpret 
the more open structures that we find for the AOY po- 
tential as a direct result of its steep attraction and the 
resulting slow equilibration of the cluster geometry. 



C. Hierarchical assembly: Superclusters 



Figure 8. Comparison between the number of particles 
TicNn^/Ni belonging to a cluster of size Uc found in simu- 
lations and the weight fraction of particles as a function of Uc 
measured in experiments (see Figures 4 and 6 of Ref. [40l|V 
The simulation results are for starting droplet size ad/uc ~ 4, 
and Tjc = 0.0034. In the experiments the amount of building 
blocks added via the water or oil phase was 108 mg. For both 
simulation and experiments the percentage of particles per oil 
was 3.4%. 



is viable. Zerrouki et al. [33| ], at 5% in weight of silica 
microparticles per volume of oil, found a fraction with re- 
spect to singlets of 0.75, 0.7 and 0.47 for doublets, triplets 
and quadruplets, respectively. Although these are higher 
yields than what we found in our simulations and exper- 
iments, their experimental parameters differ significantly 
from ours, so that no conclusion about the relative per- 
formance of both methods can be drawn. 

In simulations, we find that clusters with the same 
number of constituent particles can still have a variety of 
different structures (isomers). Instead of distinguishing 
between all possible isomers we classify clusters based on 
their number of bonds. The bond-number rif, is defined 
as the total number of bonds in a cluster and, although 
unable to distinguish between all possible isomers, gives 
an indication of the compactness of the cluster; for a 
given value of ric a smaller number of bonds indicates a 
more open structure as compared to a cluster with more 
bonds. 

Fig. [5^ shows a stacked histogram of the number of 
clusters with a specific bond-number (SWY potential, 
CTd(O) = 8ctc, rjc = 0.0034). The total height of the 
columns indicates the number of cluster N„^ with ric 
particles. Each bar is divided in differently colored re- 
gions with a relative size proportional to the number 
of clusters with ni, bonds. Each region is labeled with 
the actual bond-number. For ric = 2, 3, only one type 
of cluster is found with n^, = 1,3, respectively. Clearly 
these bond- numbers correspond to doublets and triplets, 
respectively. For ric = 4, two different structures are 
found, a small fraction of clusters with an open structure 
with only four bonds, and a structure with six bonds, 
corresponding to the tetrahedron shown in Fig. [5]3. For 



In order to investigate the possibility of hierarchi- 
cal assembly of colloidal particles using the droplet- 
evaporation technique we carried out computer simula- 
tions of a mixtures of tetrahedral clusters and emulsion 
droplets. Hence, we prepared an initial configuration 
of the simulation consisting of clusters with tetrahedral 
symmetry (Fig. [SJd). No single particles or other than 
tetrahedral cluster structures were present. 

The tetrahedral clusters are thermal in the sense that 
the particles forming the clusters are kept together solely 
by the short-ranged attraction and can in principle dis- 
solve, i.e. bonds can break on a long time scale by ther- 
mal activation. In experiments, bond breaking is even 
more unlikely because clusters are held together by van 
der Waals interactions, which are much stronger than the 
attraction used in our model. Nevertheless, structures 
that require bond-breaking in order to form are easily 
recognized. 

Figure [10] shows the structures obtained in the simu- 
lations (SWY potential, = 9(Tc, ?7c = 0.1). In partic- 
ular. Fig. llOb -b) show two structures consisting of two 
tetrahedral building blocks. The octahedral dipyramid 
(Fig. [TOb ) is formed by two tetrahedra rotated by 30° 
against each other and with two faces touching. Fig.llOb) 
shows two truncated hexagonal layers. This cluster for- 
mation is possible only because one of the initial tetra- 
hedral clusters has dissolved. This structure is therefore 
not accessible experimentally when non-thermal clusters 
are used. Fig. [TUb l shows superclusters of three tetra- 
hedra, while Fig. [TUli) shows a supertetrahedron, i.e. a 
cluster formed by four tetrahedra arranged a tetrahedral 
geometry. 



IV. SUMMARY AND CONCLUSIONS 

We investigated cluster formation via emulsion droplet 
evaporation with computer simulations and experiments. 
We used Metropolis Monte Carlo simulations to model 
the process of cluster formation in a binary mixture of 
colloidal particles and emulsion droplets. The colloidal 
particles interact via both a short-ranged attraction and 
a long-ranged repulsion, while the second component that 
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Figure 9. Number of cluster A'^ with ric colloidal particles 
(crd(O) = Sctc and rjc = 0.0034). The total height of the 
columns indicates the number of cluster N„^ with ric parti- 
cles. Each bar is divided in differently colored regions with 
a relative size proportional to the number of clusters with Ub 
bonds. The numerical label indicate the bond-number for the 
region. For clarity, the Uc — 1 bar is omitted, a) Results for 
the SWY potential. Shown are also the names of the most 
relevant structures, b) Results for the AOY potential. 



represents the emulsion droplets interacts only with an 
attractive well with the colloids. This potential well has 
a minimum at the droplet surface in order to induce 
the Pickering effect. The droplet-droplet interaction is a 
hard-core interaction with a hard-sphere diameter chosen 
so that the droplets cannot merge. The droplets shrink 
at a fixed rate, in order to model experimental conditions 



of droplet evaporation. 

We also performed experiments on polystyrene spheres 
154 nm in diameter in a toluene- water emulsion. The 
emulsion was vitrified and analyzed with cryo-FESEM, 
before the evaporation process. The micrographs indi- 
cated a random distribution of the positions of the parti- 
cles that are trapped at the droplet surface. These results 
can be surprising because charged particles act as elec- 
tric dipoles when trapped onto droplets due to the effect 
that the part of their surface exposed to the nonpolar 
solvent (toluene) cannot sustain its charge. The result- 
ing long-ranged dipolar repulsions may result in regu- 
lar spatially separated arrangements of the particles at 
the interface [5l|, [52|. For micron-sized particles this is 
supported by optical micrographs that indicate polyhe- 
dral arrangements when a small number of particles are 
bound onto an emulsion droplet [s^ . A different scenario 
can be expected for our submicron-sized particles because 
their higher diffusivity can interfere with repulsive inter- 
actions, and suppress regular orientation. Interestingly, 
an earlier study of droplets stabilized by a large number 
of submicron-sized colloids could even demonstrate that 
the particles are not necessarily kept separated from each 
other and can form close-packed islands or even a mono- 
layer [Hi. 

Computer simulation snapshots and radial distribution 
functions were used to analyze the dynamics of cluster 
formation in computer simulation. In agreement with 
experimental results we find that in our model the parti- 
cles can freely diffuse on the surface of the droplet before 
the evaporating droplets force particle agglomeration into 
clusters. The degree of ordering of the particles on the 
droplet surface depends on the range of the repulsive in- 
teraction, (Debye screening length) that in our model is 
of the order of two particle diameters. Choosing a longer 
range for the repulsion could lead to ordered distribu- 
tion of particles on the surface of the droplets. Further- 
more the strength of dipolar interactions, neglected in 
our present model, could be relevant. 

After the complete evaporation of the droplets we find 
stable clusters that range from sphere doublets to com- 
plex polyhedra. The structures and size distributions 
found in simulations matched those found in experi- 
ments. Histograms show that larger clusters can be ob- 
tained by increasing the initial size of the droplets or by 
increasing the density of colloidal particles at the expense 
of a smaller yield for smaller clusters in accordance with 
the results of Wagner et al. . 

The bond number was used to distinguish different 
structures with the same number of constituent parti- 
cles. We found that although the AOY potential gives 
the same clusters and size distributions as the SWY po- 
tential, the AOY interaction results in a larger number of 
possible structures than the SWY interaction. In partic- 
ular, the AOY potential gives structures with a smaller 
number of bonds, i.e. with more open structure. This is 
intuitively reasonable, since the steep attractive part of 
the AOY potential results in a difficult equilibration of 
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a) b) c) d) 

Figure 10. Supercluster structures (SWY potential, ad{0) — 9(7c, 77c = 0.1). Particles with the same color belong to the same 
initial tetrahedral building block, a) Octahedral dipyramid, b) truncated hexagonal layers, c) Supercluster composed of three 
tetrahedra, d) Supertetrahedron. 



the geometric structure of the clusters. Although these 
potentials do not model quantitatively our experimental 
system, this study can give an indication of what type of 
interactions one should use in order to change the cluster 
morphology to more open (softer) clusters. 

Our simple model reproduces the experimental results 
accurately despite a lack of realistic energy or time scales. 
It is therefore sensible to assume that the model captures 
the essential physics of the assembly process and that 
more complex assembly processes can be studied with a 
certain confidence. Hence, the model can be useful to 
guide experimental work. As an example we applied the 
theoretical model to a fluid mixture of tetrahedra clus- 
ters and droplets. The computer simulations show that 
the assembly process via emulsion droplet evaporation 
can lead to stable superclusters with two, three and four 
tetrahedral building blocks. These novel structures are 
not found in the assembly of single nanoparticles and 
could represent a step in the direction of novel and com- 
plex mesoscale materials. 



We neglected in our theoretical model the formation 
of dipole moments for particles trapped at the droplet 
surface. Including dipolar interactions and increasing the 
range of the repulsion both constitute interesting steps 
beyond the current work. Furthermore, using our present 
model for studying gelation, as was recently reported in 
colloidal dispersions with a small immiscible liquid (ssj . 
could be interesting. 
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